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We consider a minimal grand unified model where the dark matter arises from non-thermal decays 
of a messenger particle in the TeV range. The messenger particle compensates for the baryon 
asymmetry in the standard model and gives similar number densities to both the baryon and the 
dark matter. The non-thermal dark matter, if massive in the GeV range, could have a free-streaming 
scale in the order of 0.1 Mpc and potentially resolve the discrepancies between observations and 
the ACDM model on the small scale structure of the Universe. Moreover, a GeV scale dark matter 
naturally leads to the observed puzzling proximity of baryonic and dark matter densities. Unification 
of gauge couplings is achieved by choosing a “Higgsino” messenger. 


INTRODUCTION 

The standard model of particle physics is believed to 
be incomplete. For decades the strongest arguments are 
based more on aesthetic reasonings than on empirical 
evidence. One example is the fine-tuning in the mass of 
the scalar Higgs requires new physics at around 1 TeV to 
stabilize the electroweak scale, for which the benchmark 
solution is weak scale supersymmetry (SUSY). Another 
example is the gauge coupling unification which suggests 
a grand unified theory (GUT) at high energy scale 
0, assuming a desert between the GUT scale and the 
electroweak scale where all the particles come in complete 
multiplets of the GUT group. It turned out that in the 
minimally supersymmetric standard model (MSSM) [^, 
gauge couplings unify to a much better precision than in 
the standard model, which adds to the attractiveness of 
the weak scale SUSY. 

On the empirical side the situation has dramatically 
improved over the last few years due to insights from 
precision cosmological observations, including the exis¬ 
tence of dark matter, the acceleration of the cosmic 
expansion, the baryon-asymmetric Universe, and a nearly 
scale-invariant density fluctuations, none of which can 
be explained within the standard model. (A minimal 
model addressing all these issues has been proposed in 
1^.) Emerging from the observations is a description 
of the Universe based on the AGDM model Q: a tiny 
cosmological constant plus cold dark matter. 

Neither the MSSM nor the AGDM model is perfect, 
however. In the MSSM the most notable problems are 
the SUSY flavor problem, new flavor violations from 
various superpartners, and the non-observation of a light 
Higgs as well as any light sparticles, which implies fine- 
tuning at a few percents level Q. These problems 
prompted the proposal ||] of giving up SUSY as a 
solution to the hierarchy problem. But if gauge coupling 
unification (and the neutralino dark matter) is the main 
motivation for SUSY, there is a much simpler model, 
the standard model with “Higgsinos”, which is just as 
good 0. On the other hand, there has been evidence 


suggesting, albeit not yet conclusively, inconsistencies 
between observations and numerical simulations of the 
AGDM model on the galactic and sub-galactic scales [|j; 
it seems that ACDM model predicts too much power 
on the small scales. One way out is to introduce 
warm dark matter (WDM) |^, ^ which has a small 
free-streaming scale Aps ^ 0-1 Mpc, hence suppressing 
density fluctuations on scales smaller than Aps- 

In this letter, using unification instead of naturalness 
as the main incentive, we show that the three otherwise 
independent aspects: gauge coupling unification, the 
puzzling proximity of Udm and U{,, and the small 
free-streaming scale, can all be intertwined in a non¬ 
trivial way. It is based on the model discussed in |l^ 
where the dark matter arises from late-time decays of 
a heavy messenger particle compensating for the baryon 
asymmetry in the standard model. First, gauge coupling 
unification suggests the existence of “Higgsinos,” which 
is the messenger particle, at the TeV scale and sets 
the GUT scale to be ~ 10^"* GeV [Q. Next cosmology 
constrains the decay temperature of the “Higgsinos” to 
be ^ 10 MeV and determines the scale of the higher 
dimensional operators responsible for the Higgsino decay 
to be GeV |l^, a most natural value in the GUT 

picture. Then the measured ratio of HDivi/Hb can be 
used to fix the dark matter mass to be ~ I GeV, which 
turns out to come with the necessary free-streaming 
scale to suppress small scale structures of the Universe. 
Alternatively, the need for a sub Mpc free-streaming scale 
can be used to argue for a GeV scale dark matter, with 
which the Hdm is naturally close to the fib- 


THE MODEL 

The basic idea in |ll| is that the dark matter S is 
produced non-thermally from the decay of a heavy mes¬ 
senger particle X, which carries the baryon number and 
compensates for the baryon asymmetry in the Universe. 
Both S and X, which we call the dark sector, are assumed 
to be charged under a Z 2 symmetry, the T-parity, while 
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the whole standard model is T-even. The dark matter S 
is then a stable particle being the lightest T-odd particle 
(LTP). At the time of baryogenesis, we assume that the 
B — L number is distributed between the T-even and 
T-odd sectors, resulting in the following relation: 

= -qB-L{nx - n^) , (1) 

where qs-L is the B — L charge of the messenger A, and 
and n^_j^ are the B — L number densities in the 
standard model and the dark sector, respectively. On 
the other hand, since both X and X eventually decay 
into the LTP, the dark matter candidate S, its number 
density is given by the total number of X and X particles 


n-DM — = nx + tix, (2) 

which is independent of the in Eq. (|^) and would 

suggest there is no connection between the baryonic and 
dark matter densities, unless nx ^ nx 0 or the other 
way around. This implies the lifetime of X should be long 
enough so that it does not decay until after most of the X 
particles annihilate with A, which will be the case if there 
is no relevant or marginal operator contributing to the 
decay of A. Then at temperature T < mx, where mx is 
the mass of the messenger, particle A starts to annihilate 
with its anti-particle A through gauge interactions and 
we are left with an abundance of A. Consequently, 

~ ~'^B-L — ~QB-L ■'T'DM • (3) 

This is a very general framework. In | [TI| we worked 
out the cosmological constraints, as well as an example 
and the associated collider phenomenology. The model 
we are interested here, is that with a singlet scalar dark 
matter S and a “Higgsino” messenger A. Recently it 
is pointed out 0 that standard model plus Higgsinos 
achieves gauge coupling unification at around 10^^ GeV 
with an accuracy similar to that of the MSSM. At two- 
loop level, for mx = 1 TeV this minimal model predicts 
smaller values of Ug than measured by 0(5%), whereas 
larger values are predicted in the MSSM. A GUT scale 
in the order of 10^^ GeV is too small for a conventional 
GUT model due to the constraint from proton decay. 
One possibility, as discussed in |Q , is to embed the setup 
in extra-dimensional or deconstructed models to lower 
the string scale down to 10^“^ GeV or so. This has the 
bonus of resolving the triplet-doublet splitting problem 
in conventional GUT models. A crucial departure from 
the model in |Q arises in that there the dark matter is a 
mixture of the neutral component of the Higgsinos and an 
extra singlet fermion, which is the usual WIMP scenario, 
whereas in our model the dark matter comes from non- 
thermal decays of Higgsinos and has a mass which will 
be determined to be in the GeV range. 

The scenario proceeds in three stages. In the first 
stage, baryogenesis is achieved by the out-of-equilibrium. 


GP-violating decay of a T-odd scalar particle P into 
ii X and -|- A, where ti is the lepton doublet in 
the Ah generation in the standard model. CP violation 
enters through the phases in the Yukawa couplings of 
the P with ii and A if there are more than one T’s. So 
at one loop one obtains an asymmetry proportional to 
the imaginary parts of the Yukawa couplings of the P’s. 
After the P’s drop out of the thermal equilibrium, we are 
left with an asymmetry in lepton number in the standard 
model. Note that this mechanism is reminiscent of the 
leptogenesis |Q, in which T is a heavy neutrino and A 
is the ordinary scalar Higgs. The effective B — L number 
of A in Eq. (||) is determined to be qB-L = —1- 

The second stage is the complete annihilation of the 
Higgsino with its anti-particles. This constrains the mass 
of the messenger particle mx so that nx ^ nx ~ 0. 
The thermally averaged cross section is estimated to be 
/{2i)inTm?x) for s-wave annihilations into SU{2) gauge 
bosons which gives an upper bound on mx- 

mx < 150 TeV (. (4) 
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However, in order to establish the connection between 
the baryon and dark matter number densities, n^^j^ ^ 
—qB-LnjoM, nnx is preferred to be much lower than the 
upper bound |n . Moreover, a heavy A spoils the gauge 
coupling unification. We found that, increasing mx from 
1 TeV to 100 TeV, the accuracy is worsened by 0(1%). 
On the other hand, the lower bounds on mx simply 
come from direct searches of new fermions which is about 
100 GeV. Therefore it is quite natural for the messenger 
particle to have a mass in the 0(TeV) range. 

The third stage is the late time decay of the Higgsinos 
into the dark matter. The decay has to occur before 
the big bang nucleosynthesis (BBN), and after the com¬ 
pletion of the pair annihilation, as previously explained, 
which requires the Higgsino to survive long after it falls 
out of thermal equilibrium: 


10 MeV < T, « ^ , (5) 

where T^j is the temperature at which A decays. The 
long lifetime can be realized by assuming that A only 
decays through a dimension five operator 

0decay = ^(V5)(ife^), (6) 

where H is the scalar Higgs and e° is the right-handed 
charged lepton in the Ah generation. The above bounds 
on Td translate into bounds on M: 


10^° 



<c M < lO^"* 


GeV 



3/2 


( 7 ) 

If we assume that the particle responsible for baryogene¬ 
sis, T, is thermally produced, the reheating temperature 
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of the Universe is restricted to be 

™x<mp<T«<lTev(^^^^) , (8) 

where the upper bound comes from the requirement 
that the thermally produced S through the operator in 
Eq. is less significant than the non-thermal com¬ 
ponent. Eq. (H) gives a stronger lower bound on M: 
10^^ GeV (mx/lTeV)^/^. Therefore, to be consistent 
with Eq. we need mx ^ 1 TeV. For example, if 
mx 10 TeV, M can be ~ 3 x GeV and 

the reheating temperature could be 10 — 10^ TeV. We 
emphasize that the lower bound in Eq. (|^) depends on 
the assumption that the P particle is thermally produced. 
In a scenario where P dominates the energy density of 
the Universe and then decays into X and radiation, Tr 
can be much smaller than mx- Then mx is bounded 
below only by direct searches. 

From the discussions above, we see that it is quite 
natural to identify M with Mqut ^ 10^"^ GeV. In this 
case, the decay temperature Td is just before the BBN, 
Td ~ 10 MeV, for mx ~ 1 TeV. A priori it is not 
clear at all that the cosmological bounds on M should be 
consistent with the GUT scale ~ 10^^ GeV in our model. 
For example, if the messenger were to decay through a 
dimension six operator, then the bounds on M would lie 
between 10® and 10^ GeV 0 , which would make the 
identification with Mgut less convincing. 

There are in fact two dangerous marginal operators 
allowed by T-parity: O 4 = liXS, which contributes to 
the Higgsino decay, and O 4 = SSH^H, which could 
thermally produce S. The first one, being a Yukawa- 
type coupling, will be absent if we set the initial value 
at Mgut to zero, which we will do. The second one 
can be fine-tuned away, since we are not motivated by 
naturalness principle. In fact, if we extend the model 
a little bit by using a complex scalar for 5, then O 4 
is forbidden by a Z^, symmetry under which both S 
and X have unit charge. Then the Higgsinos decay 
through X£SS. (Another dangerous operator, X£S\ 
can be forbidden by a Z 2 symmetry under which X is 
even and S is odd.) In addition, O 4 = SS'^HH^ is 
generated radiatively with a coefficient suppressed by 
lepton Yukawa couplings, thus alleviating the fine-tuning. 
In any case, the scalar masses, the Higgs and the dark 
matter, are fine-tuned. 

Note that even though X decays into dark matter 
plus leptons, the excess in lepton number doesn’t get 
converted into the baryon number because the decay 
always happens after the electroweak phase transition 
[0 and the sphaleron process is ineffective. We then 
obtain the number densities ns and tidm as follows: 


UB = e , 


noM = 


(9) 


where the efficiency e is the relation between the B — L 
number and the baryon asymmetry in the presence of the 


sphaleron process. It is different from the standard model 
value , now that the presence of X modifies the charge 
neutrality conditions, and calculated to be 25/79 with 
the additional constraint that the total {B — L) number 
is zero before the decay of X. Therefore the ratio of the 
baryonic and dark matter densities UDivi/Hf, is given by 


Hdm 

rib 


3.16 


todm 

TOp 


( 10 ) 


from Eq. i) , where mp is the proton mass. The measured 
ratio is SlDM/Hf, ~ 5.1 which implies muM 1.6 rup = 
1.5 GeV. Eq. ( |l^ is a central prediction of our model: the 
baryonic mass density in our Universe will be close to the 
dark matter density if the dark matter has a mass close to 
the nucleon mass. Indeed, we will see that considerations 
from reducing the small scale structure of the Universe 
prefers light dark matter, C>(GeV), showing a non-trivial 
consistency in this framework. 


COSMOLOGY 

The cosmology with a dark matter candidate arising 
from non-thermal decays is very interesting. In the 
ACDM model the dark matter decouples from the ther¬ 
mal bath while being non-relativistic and has a velocity 
distribution peaked at around zero at the time of struc¬ 
ture formation. In other words, the momentum profile is 
Maxwellian. On the other hand, the non-thermal dark 
matter (NTDM) is relativistic when produced through 
decays of the messenger particle and its momentum 
distribution is peaked at tox/ 2, which subsequently gets 
red-shifted by the expanding Universe and becomes non- 
relativistic. 

Such a primordial velocity dispersion has important 
implications in structure formation, as the dark matter 
can smooth out inhomogeneities by streaming out of 
overdense regions and into underdense regions, the so- 
called Landau damping |Q, which is characterized by 
the free-streaming scale Aps- By definition cold dark 
matter has a small Aps that is irrelevant for structure 
formation, whereas hot dark matter, Aps ^ 40 Mpc, 
has such a large free-streaming scale that would prevent 
galaxy formation in the early epoch. So far observations 
on large scale structures clearly prefers a dark matter 
candidate that is cold, and hence the AGDM model. 
Nevertheless, as mentioned, on small scales (sub Mpc) 
simulations of GDM seems to be at odds with observa¬ 
tions. First, the substructure of GDM halos is predicted 
to be richer than observed. Second, the simulated 
density profiles of dark matter halos are generally cuspier 
than inferred from rotation curves. At this moment 
it is not clear if this is indeed a failure of the ACDM 
model, as there are complex issues with large N-body 
and hydrodynamics simulations. 
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These discrepancies do not diminish the tremendous 
successes of the ACDM model on larger scales, however, if 
these problems are real, they present a great opportunity. 
The simplest known mechanism for smoothing out small 
scale structures is the Landau damping. In this regard a 
popular candidate is the WMD, hot dark matter cooled 


down, which has a free-streaming scale |10 


Afs = 0.2 


2. 1 /mwDM\ 3 


V keV / 


Mpc. (11) 


Observed properties of Lyman a forest constrain the 
power spectrum at small scales and therefore put a lower 
bound on the mass of the WMD ||^: towdm > 750 eV, 
which translates into Afs < 0.16 Mpc for OwDMh^ = 
0.15. (A more recent analysis puts a weaker bound: 
mwDM > 550 eV [^.) Furthermore, high-resolution 
cosmological N-body simulations seem to find better 
agreements with observations for Afs ~ 0.1 Mpc 0 - 
For comparison, we can calculate the free-streaming scale 
for the NTDM in our case [l^ . 



where v{t) and a{t) are the physical velocity and the 
FRW scale factor, respectively, of the NTDM at time t, 
and ^eq denotes the time for matter-radiation equality. 
Therefore, for a prototypical mx/nis/Td = 3 TeV/1 
GeV/10 MeV scenario that is well-motivated from previ¬ 
ous discussions, the NTDM is able to produce enough 
power on the small scales to be consistent with the 
Lyman a forest data and potentially resolve the dis¬ 
crepancies mentioned above. There are also constraints 
coming from studies of phase space density Q, which 
are weaker than those coming from Lyman a forest . 

On the other hand, the NTDM is not exactly the same 
as the WDM, since they still have different momentum 
distributions. As an example, power spectrum of non- 
thermal production of neutralinos by the decay of topo¬ 
logical defects was considered in jl^ and found to be 
different from that of a 1 keV WDM for k > 5 /iMpc”^. 


SUMMARY 

Using gauge coupling unification and cosmology as the 
main hints for physics beyond the standard model, we 
have considered a minimal GUT model which attributes 
a common origin to the baryon asymmetry and the 
dark matter, giving similar number densities to both 
the dark matter and the baryon. The dark matter, 
a singlet scalar, is produced by non-thermal decays of 
a messenger particle, the Higgsinos, whose existence 
implies the unification at ~ 10^"* GeV. There are checks 


on the the model from several orthogonal directions. 
Gosmological bounds point to a Higgsino in the TeV 
range and a mass scale consistent with the GUT. The 
ratio DDivi/f^b can be used to determine the mass of the 
dark matter to be GeV, which in turn gives a sub Mpc 
free-streaming scale consistent with observations and, at 
the same time, has the potential of resolving the ACDM 
crisis by reducing the power spectrum on small scales. 
Cosmology may be the best arena to test this model. 
The Higgsinos, on the other hand, might be too heavy to 
be discovered in the coming collider experiments. 

Acknowledgments We benefited from conversations 
with Neal Dalai and Carlos Peha-Garay. This work is 
supported by funds from the IAS and in part by the 
DOE grant number DE-FG02-90ER40542. 


[1] H. Georgi and S. L. Glashow, Phys. Rev. Lett. 32, 438 
(1974); H. Georgi, H. R. Quinn and S. Weinberg, Phys. 
Rev. Lett. 33, 451 (1974). 

[2] S. Dimopoulos and H. Georgi, Nucl. Phys. B 193, 150 
(1981); S. Dimopoulos, S. Raby and F. Wilczek, Phys. 
Rev. D 24, 1681 (1981). 

[3] H. Davoudiasl, R. Kitano, T. Li and H. Murayama, 
arXiv:hep-ph/0405097. 

[4] See, for example, D. N. Spergel et al. Astrophys. J. Suppl. 
148, 175 (2003) [arXiv:astro-ph/0302209]. 

[5] L. Giusti, A. Romanino and A. Strumia, Nucl. Phys. B 
550, 3 (1999) [arXiv:hep-ph/9811386]. 

[6] N. Arkani-Hamed and S. Dimopoulos, arXiv:hep- 
th/0405159. 

[7] N. Arkani-Hamed, S. Dimopoulos and S. Kachru, 
arXiv:hep-th/0501082. 

[8] For a review, see J. P. ©striker and P. J. Steinhardt, 
Science 300 (2003) 1909 [arXiv:astro-ph/0306402]. 

[9] S. Colombi, S. Dodelson and L. M. Widrow, Astrophys. 
J. 458, 1 (1996) [arXiv:astro-ph/9505029]. 

[10] P. Bode, J. P. Ostriker and N. Turok, Astrophys. J. 556, 
93 (2001) [arXiv:astro-ph/0010389]. 

[11] R. Kitano and I. Low, Phys. Rev. D 71, 023510 (2005) 
[arXiv:hep-ph /0411133]. 

[12] M. Fukugita and T. Yanagida, Phys. Lett. B 174, 45 
(1986). 

[13] J. A. Harvey and M. S. Turner, Phys. Rev. D 42, 3344 
(1990). 

[14] E. W. Kolb and M. S. Turner, “The Early Universe,” 
Redwood City, USA: Addison-Wesley (1990). 

[15] V. K. Narayanan, D. N. Spergel, R. Dave and C. P. Ma, 
arXiv:astro-ph/0005095. 

[16] M. Viel, J. Lesgourgues, M. G. Haehnelt, S. Matarrese 
and A. Riotto, arXiv:astro-ph/0501562. 

[17] P. Colin, V. Avila-Reese and O. Valenzuela, arXiv:astro- 
ph/0009317; V. Avila-Reese, P. Colin, O. Valenzuela, 
E. D’Onghia and C. Eirmani, arXiv:astro-ph/0010525. 

[18] C. J. Hogan and J. J. Dalcanton, Phys. Rev. D 62, 063511 
(2000) [arXiv:astro-ph/0002330]. 

[19] W. B. Lin, D. H. Huang, X. Zhang and R. H. Bran- 
denberger, Phys. Rev. Lett. 86, 954 (2001) [arXiv:astro- 
ph/0009003]. 








